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ABSTRACT 

We present 3-dimensional SPH simulations of supernova explosions from 100 
seconds to 1 year after core-bounce. By extending our modelling efforts to a 
3-dimensional hydrodynamics treatment, we are able to investigate the effects 
of explosion asymmetries on mixing and 7-ray line emergence in supernovae. 
A series of initial explosion conditions are implemented, including jet-like and 
equatorial asymmetries of varying degree. For comparison, symmetric explosion 
models are also calculated. A series of time slices from the explosion evolution 
are further analyzed using a 3-dimensional Monte Carlo 7-ray transport code. 
The emergent hard X- and 7-ray spectra are calculated as a function of both 
viewing angle and time, including trends in the 7-ray line profiles. We find 
significant differences in the velocity distribution of radioactive nickel between 
the symmetric and asymmetric explosion models. The effects of this spatial 
distribution change are reflected in the overall high energy spectrum, as well as 
in the individual 7-ray line profiles. 



Subject headings: black hole physics — supernovae: general — stars: neutron 
gamma rays: theory 



1. Introduction 

Among the many surprises that supernova (SN) 1987A brought astronomers was 
the early emission of X-rays and 7-rays (X-rays: e.g. Dotani et al. 1987; Sunyaev et al. 
1987; 7-rays: e.g. Cook et al. 1988; Mahoney et al. 1988; Matz et al. 1988). This high 
energy emission, arising from the decay of ^^Co, appeared nearly 6 months earlier than was 
predicted by theoretical models (e.g. Pinto & Woosley 1988a, Arnett et al. 1989) and led 
theorists to conclude that the ^^Ni, produced deep in the core of this exploding star, had 
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mixed into the outer layers of the supernova ejecta (e.g. Pinto & Woosley 1988b, Arnett et 
al. 1989 and references therein). Additional evidence all seems to support the mixing of 
SN 1987A's ejecta: models of the supernova light curves (Arnett 1988; Shigeyama, Nomoto, 
& Hashimoto 1988; Woosley 1988; Shigeyama & Nomoto 1990) and explanations of spectral 
line widths (Haas et al. 1990; Spyromilio, Meikle, & Allen 1990; Tueller et al. 1990). 
SN 1987A is not peculiar in this mixing. Many supernovae show evidence of mixing in their 
spectra (e.g. Spyromilio 1994, Fassia et al. 1998) and the hght curves and spectra of Type 
lb SN seem to be best fit by mixed models (Shigeyama et al. 1990; Woosley & Eastman 
1997). It appears that mixing is a generic process in core-collapse supernovae. 

These results have stimulated a series of multi-dimensional hydrodynamical simulations 
trying to produce the observed mixing (Arnett, Fryxell, & Miiller 1989; Hachisu et al. 
1990, Pryxell, Arnett & Miiller 1991; Miiller, Pryxell & Arnett 1991; Herant & Benz 
1991; Herant & Benz 1992; Herant & Woosley 1994; Kifonidis et al. 2000; Kifonidis et 
al. 2003). Although these simulations seem to be able to explain the mixing in Type lb 
supernovae (Kifonidis et al. 2000), none of these results are able to explain the extensive 
mixing observed in SN 1987A. Possible solutions to this shortcoming have been proposed: 
(1) perhaps the decay of ^^Ni injects enough energy to force additional mixing (Herant & 
Benz 1992) or (2) convection in the pre-coUapse core provides enough seeds to enhance 
mixing (Herant & Benz 1992). A third possibility is that the supernova explosion itself 
is asymmetric (Nagataki et al. 1997; Nagataki et al. 1998; Nagataki 2000 and references 
therein). Nagataki et al. (1998) found that not only could shght asymmetries in the 
supernova explosion produce the required mixing to explain 1987A, but they could also 
explain anomalies in the nucleosynthetic yields produced by several supernovae. 

The evidence for such global asymmetries in supernova explosions has been growing 
steadily. One piece of evidence arises from attempts to understand the high space velocities 
of neutron stars. The high observed velocities of pulsars, along with evidence of neutron 
star/remnant associations, and the formation scenarios of neutron star binaries all suggest 
that neutron stars are given strong kicks at birth. These kicks are most easily explained 
by some asymmetry in the supernova explosion (see Fryer, Burrows, & Benz 1996 for a 
review). In addition, the most straightforward explanation of the large polarization seen in 
core-collapse supernovae (see Wang et al. 2001; Leonard & Filippenko 2001 and references 
therein) is that the explosion driving these supernovae is inherently asymmetric (Hoflich 
1991). 

The asymmetries are believed to have their origin in the explosion mechanism itself. 
For instance, even if the collapsing star is initially spherically symmetric, some asymmetry 
can be produced due to convection taking place in, and above, the proto-neutron star 
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(Herant et al. 1994; Burrows, Hayes, & Pryxell 1995; Janka & Miiller 1996). To date, 
these asymmetries are not extreme enough in the theoretical models to explain the mixing. 
Large asymmetries may occur if the collapsing star is asymmetric due to nuclear burning 
(Burrows & Hayes 1996; Lai 2000), however, sufficiently large departures from spherical 
symmetry have only been produced by assuming extremely asymmetric collapsing cores 
(Burrows & Hayes 1996). Alternatively, rotation can produce strong asymmetries in the 
supernova explosion (Monchmeyer & Miiller 1989; Janka & Monchmeyer 1989; Fryer & 
Heger 2000, Khokhlov et al. 1999). The nature of these asymmetries depends upon the 
angular momentum profile of the collapsing star and, although most calculations predict 
jet-like explosions along the rotation axis, some calculations imply that an equatorial 
explosion could occur (Monchmeyer & Miiller 1989). 

In this paper, we follow supernova explosions in 3-dimensions from 100 seconds to 1 
year after core-bounce. We model a series of initial explosion conditions with both jet-like 
"axial" and equatorial asymmetries of varying degree. In §2, we describe these simulations 
and their results with comparisons to past work. In §3, we discuss the 7-ray emission from 
these explosions and present calculations of 7-ray spectra as a function of time and viewing 
angle. We conclude with a disscussion of the observational prospects of these results, 
highlighting upcoming gamma-ray missions. 

2. Explosion Simulations 

For our hydrodynamic simulations, we have used the 15 Mgprogenitor (sl5s7b) by 
Weaver & Woosley (1993). This star has been evolved with a piston-driven explosion to 
100 s after bounce, producing 0.24 Mgof ^^Ni. The total energy of this model is roughly 
1.5 X 10^^ erg with roughly 1.0 x 10^^ erg in kinetic energy. As this explosion moves through 
the star, the shock hits composition boundaries where strong entropy gradients exist. When 
the shock hits these boundaries, Rayleigh- Taylor instabihties develop, which can grow and 
cause the star to mix (Chevaher & Klein 1978; Weaver & Woosley 1980). Our simulations 
model this mixing and concentrate on the effects that asymmetries have on it. 

2.1. Numerical Schemes 

We model a series of explosions from 100 s to lyear after the launch of the supernova 
shock (Table 1) using a 3-dimensional smooth particle hydrodynamics code (see Warren et 
al. 2002 for details) based on the parallel oct-tree algorithm developed by Warren & Salmon 
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(1993). This lagrangian code tracks the composition exactly. For 2 of our simulations, 
we included the energy injection from ^^Ni and ^^Co decay (and hence also trace the 
abundances of ^^Co and ^^Fe). For these 2 models, we assume that all of the decay energy 
is injected into the gas. At late times, this will overestimate the total energy deposited, 
as a fraction of this energy will escape, but since we would like an upper limit on the 
effects of ^^Ni decay, and because the injection of energy is less important to the explosion 
dynamics at late times, this assumption is adequate. For decay energy, we assume that 
the total energy from decay is 9.3 x lO^^ergg"^ with 33% being released from ^^Ni decay 
with a 6. Id half-life and the other 67% being released from ^^Co decay with a 77d half-hfe 
(Colgate, Petschek, & Kriese 1980). This is similar to what Herant & Benz (1992) and 
Herant & Wooslcy (1994) used for their models. For an equation of state, we use the same 
"low-density" equation of state developed by Blinnikov, Dunina-Barkovskaya, & Nadyozhin 
(1996) that is used in our progenitor model. 

We mapped the Weaver & Woosley (1993) model with its 100 s long, spherically 
symmetric explosion into our 3D SPH code with 2.2 million variably-massed particles. We 
model the entire sphere, so there are no axis boundaries in this simulation and the neutron 
star mass at the center is added through an external gravitational force. Asymmetries 
are added artificially to the velocities only and, for most of the explosions, are added in 
such a way to insure that the total kinetic energy of the explosion is conserved. The two 
asymmetries we model are jet explosions: 

^^radial =(« + /? X \z\/rXl^^^^ (1) 

and equatorial explosions: 

^^radial = (« - /3 X \x\ / r)vl^J^^^ (2) 

where v^^^^^^ is the velocity from the one dimensional calculation, Wradiai is the radial velocity 
for the asymmetric setup, and x, z, r are the x-position, z-position, and radius of the 
particle. The values of a and /3 for each model are given in Table 1, along with the 
initial kinetic energy of each model. The magnitude of the asymmetries are guided by the 
results of Fryer & Heger (2000), who found that their rotating core-collapse simulations in 
2-dimcnsions had velocities that were a factor of 2 higher in the pole than in the equator 
1.5 s after bounce. We do not use the extreme asymmetric explosions of Khokhlov et al. 
(1999) which concentrate the explosion energy into a narrow jet. As we shall show (and 
as Nagataki 2000 argued), such extreme asymmetries are not required to explain outward 
mixing of nickel in supernovae. Although at the end of their simulation, the asymmetry in 
polar vs. equatorial velocities is growing, the trend as the shock moves outward will be to 
spherize the shock. So at this point, it is not clear whether the shock at 100 s will be more 
or less asymmetric than what Fryer & Heger (2000) found at the end of their simulation. 



- 5 - 



After mapping these models into our 3-dimensional SPH code, we then run the 
explosion out to 1 year. When the shock reaches the edge of the star, we assume it is moving 
through a vacuum. In reality, such a star will have a stellar wind atmosphere surrounding 
it, but the density of a ISMgWind is so low, that for the purposes of our simulation, 
zero density material is appropriate. In addition, photon transport is not modeled in our 
calculations. However, prior to shock breakout, the photons are essentially trapped in 
the shock. By the time the shock breaks out, most of the internal energy has already 
been converted to kinetic energy, so although including photon transport will change 
the ionization state of the exploding stellar material, it does not affect the kinematics 
significantly. 

The smooth particle hydrodynamic technique automatically adds a level of perturbation 
in the code. With our initial conditions, these perturbations are random with a maximum 
1 — (7 deviation in a shell of 5-7%. Such large deviations only occur at composition 
boundaries, specifically the silicon and oxygen shells, and match well the deviations arising 
from explosive oxygen and sihcon flashes prior to collapse (e.g. Bazan & Arnett 1998). 

2.2. Explosion and Nickel Distribution 

Although we use the same 15 Mgprogenitor that was used by Herant & Woosley 
(1994), it has a higher explosion energy and we map this model onto our 3-dimensional grid 
200 s earlier than Herant & Woosley (1994) mapped their spherically-symmetric explosion 
simulation onto a 2-dimensional grid. Therefore, although their study has the closest 
similarities with our work, it is difficult to make direct comparisons to their simulations. 
Nevertheless, it is interesting to compare the velocity distribution of each chemical element 
from our 3-dimensional simulations with the 2-dimensional simulations of Herant & Woosley 
(1994) at 90 days (compare the lower right panel in Fig. 1 of Herant & Woosley 1994 with 
Fig. 1 in this paper). In Fig.l, the material labelled "hydrogen" includes all material in the 
hydrogen envelope (as did Herant & Woosley 1994). Similarly, by "nickel" we refer to both 
the distribution of nickel as well as its decay products (most notably ^^Co) . Although the 
distribution of elements is similar in both the 2 and 3 dimensional simulations, the stronger 
3-dimensional explosion causes all of the ejecta to be moving slightly faster than that of the 
2-dimensional simulation and it is difficult to compare mixing instabilities. 

Comparing the convective instabilities themselves is also difficult. In 3-dimensions, 
the "mushroom" -like structures formed by Rayleigh- Taylor instabilities are not so well 
defined, and don't lie along any one plane. However, 4.3 hours into the explosion, it is 
clear that instabilities have developed (Fig. 2) and these instabilities ultimately mix nickel 



-6- 



knots far out into the star (Fig. 3). At the start of the explosion, ^^Ni is found only 
in the inner l.GMgof the star (1.3 Mghecomes the neutron star, so the nickel is limited 
to the inner 0.3MQof ejecta). By the end of the simulation, ^^Ni has mixed out nearly 
to 5Mq, beyond the ~4.5 Mghoundary that marked the edge of the helium layer (Fig. 
4). Unfortunately, this mixing is, if anything, less than the amount of mixing found in 
the 2-dimensional simulations of Herant & Woosley (1994). The fact that the mixing is 
less in 3-dimensions vs. 2-dimensions could be due to the lower effective resolution (we 
only have 2.2 million particles in 3-dimensions vs. 25,000 particles in the 2-dimensional 
simulations)^. However, bear in mind that the turbulent inverse cascade behaves differently 
in 2 and 3 dimensions (2-dimensional inverse cascades drive energy to large scales whereas 
3-dimensional simulations suggest the energy is driven to small scales and dissipated). It is 
likely that these differences cause the 2-dimensional simulations to produce more extended 
instabilities. In any event, it appears that neither our spherical 3-dimensional simulations 
nor the 2-dimcnsional simulations seem to give enough mixing to explain the observations 
of supcrnovac like 1987A. 

Nagataki et al. (1998) and Nagataki (2000) found that they required mild asymmetries 
("^-'Poic/'^Equator = 2) to explain SN 1987A. Since we model a ISMgStar, not a SN 1987A 
progenitor, it is difficult to both compare with this past work as well as constrain our results 
with observations of SN 1987A. But we can discuss the basic trends caused by asymmetries. 
Fig. 5 shows model Jet2 1 year after explosion. Note that although the density distribution 
has spherized as the shock propogates through the shallow density gradients of the red 
supergiant envelope (Chevaher & Soker 1989), the distribution of ^^Co (the decay product 
of ^^Ni) retains a large asymmetry. We discuss the effects of these asymmetries on the 7 ray 
emission in §3. 

Like Nagataki et al. (1998) and Nagataki (2000), we find that the asymmetries 
broaden the velocity profile of ^^Ni (Fig. 6). However, in our simulations, mild asymmetries 
(^Poic/^^Equator = 2) led to ouly a small increase in th maximum nickel velocity from 
2500 km s^^ to 2900 km s^^. For Nagataki et al. (1998), such mild asymmetries increase 
the maximum nickel velocity from 2200 km s^^ to 3200 km s^^! This difference could be 
progenitor dependent, an effect of 3-dimensional vs. 2-dimensional convection, or the lack of 
resolution in our 3-dimensional models. Extracting the true cause of this difference awaits 
future calculations with similar initial conditions. 

However, increasing the amount of asymmetry by another factor of 2 (Models Jet4, 



^Note that the 2-dimensional simulations of Herant & Woosley (1994) have poor mass resolution and low 
resolution could be a problem in the 2-dimensional simulations as well as the 3-dimensional simulations 
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Eq4) causes some nickel to be ejected at very high velocities. The amount of mixing in 
these cases reaches extremes with the nickel well into the hydrogen envelope. Such mixing 
has decided signatures in both the emergence of the 7-ray line fiux and the shape of the 
7-ray lines (§3). Note that the energy released from the decay of nickel also helps to mix 
out the nickel (on par with the effects of mild asymmetries). Clearly, the energy released 
from the decay of nickel cannot be neglected in any accurate mixing calculation. 

This mixing also has important repercussions for nucleosynthctic yields and the 
mass-cut for the remnant mass. Most black holes are formed in stars which produce 
supernova explosions that arc too weak to throw off all of the stellar envelope and the 
subsequent fallback produces a black hole (Fryer & Kalogera 2001). In our models (Jet2, 
Sym+Decay, etc.), more than 10% of the nickel produced is ejected well beyond the helium 
core (Fig. 7). If this trend holds for more massive stars such as the progenitor of SN 1997D 
(Turatto et al. 1998), then the entire helium core of such a star (>8Mg) could fall back 
and still enough nickel would escape to power the observed light curve! Bear in mind, 
however, that weaker explosions may well produce less mixing, so adapting the results of 
our simulations to supernovae like 1997D must be taken with some caution. 
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Table 1. Explosion Simulations 



Model^ 


^^ole / Equator 


a, (3^ 


g Initial 

(lO^^ergs) 


g Initial 

(lO^^ergs) 


K.E. 

(>Me) 


Mixing^ 


Sym 


1.0 


1,0 


1.0 


0.51 


1.3 


4.6 


Sym D 


1.0 


1,0 


1.0 


0.51 


1.3 


5.5 


Jet2 


2.0 


7377,^377 


1.0 


0.51 


1.3 


5.6 


Jet2 D 


2.0 


a/377, a/377 


1.0 


0.51 


1.3 


5.7 


Jet4 


4.0 


a/i77, a/977 


1.0 


0.51 


1.3 


8.0 


Eq2 


0.5 


4/3,2/3 


1.0 


0.51 


1.3 


5.0 


Eq4 


0.25 


8/5,6/5 


1.2^ 


0.51 


1.6 


11.3 



^The models include ones with symmetric initial conditions (Sym), explosions along the 
polar axis (Jet), and explosions along the equator (Eq). Those models which include the 
effects of Nickel and Cobalt decay are denoted with a "D" suffix. 

'^For the polar explosions, the radial velocity (fradiai) is given by: Vradiai = (« + /5 x 
l'^l/^)'"mdiai ^"^^ equatorial explosions: -Uradiai = (« - /3 x \x\ /r)v'^^^^^^-^^ where -w^Idiai ^^'^ 
velocity from the one dimensional calculation, x, z, r are the x-position, z-position, and 
radius of the particle. 

'^We limit the extent of mixing by the furthest position (in mass coordinates) beyond 
which more than 3% of the nickel produced near the core is mixed out. 

"^Note that the initial energy for the extreme asymmetric model was larger than the 
other models. In part, this explains the extended mixing of this model. 
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3. High Energy Spectral Calculations 

For our spectral calculations, we have used data from the 3-dimensional explosion 
simulation discussed in the previous section. We input ejecta material properties from five 
different snapshots in time at 150, 200, 250, 300 and 365 days after explosion. Spectral 
calculations were carried out for both the Jet2 and Symmetric explosion models. Our 
analysis of these model spectra concentrates on the differences in total luminosity and line 
profile shape with the introduction of realistic explosion asymmetries. Since the progenitor 
star used as input to our simulations was a 15 Mgred supergiant, we are unable to directly 
compare our calculated spectra with the observed high energy spectra of SN 1987A. 
However, we discuss how our models compare to various spectral trends observed from 
SN 1987A. 



3.1. Numerical Schemes 

We used a Monte Carlo technique, similar to that described in Ambwani & Sutherland 
1988, for modelling 7-ray transport in 3-dimensions. Input models of the supernova ejecta 
(element abundances, density and velocities) were taken from the "Jet2" and "Symmetric" 
SPH explosion simulations and mapped onto a 140 xl40 xl40 cartesian grid. Escaping 
photons were tallied into 250 coarse energy bins, with finer binning at the decay line 
energies to provide line profile information. The emergent photons were also tallied into 
11 angular bins (A^ = 10°) along the polar axis (the models investigated in this work are 
axisymmetric, alleviating the need to tally in azimuthal angle as well.) 

The decay of the radioactive species (predominantly ^^Ni and its decay product: 
^^Co) in the supernova ejecta gives rise to the 7-ray line emission. As in Ambwani & 
Sutherland 1988, we assign the energy of the emitted photon packets according to the decay 
probabilities (Lederer & Shirley 1978) for the various radioactive species (^^Ni, ^^Co, ^^Co, 
44rp-^ ^Sc, and ^^Na). We include a total of 56 decay lines from these species, but for the 
explosion times considered, the packets fall predominantly into ~ 15 decay lines. Roughly 
10^ Monte Carlo photon bundles were generated for each input model in proportion to the 
mass of radioactive material distributed throughout the ejecta. The material properties of 
the ejecta were not evolved with photon fiight time. However, we found that 99 % of the 
photons contributing to the observed model spectra have escape times of less than 2 days. 
This is sufficiently shorter than the timescales for change in the hydrodynamic models that 
our assumption of a fixed material background should be valid for the time slices considered 
here. The luminosity weight of each photon packet and the opacities seen by the packet 
were calculated in the comoving frame of the fiuid, but all photon properties were boosted 
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to the observer's frame before being tallied into spectral observables. 

In all models, photoelectric and pair production opacities were calculated for the 
elements H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Ti, Cr, Fe, Co, and Ni which correspond to 
the elements used in the nucleosynthesis calculations for the progenitor star from Weaver 
& Woosley (1993). The cross section data for these elements were taken from the LLNL 
Evaluated Nuclear Data Library (PIcchaty, Cullcn, and Howerton 1981, revised 1987). The 
angle- and energy-dependent Compton scattering opacities were calculated assuming that 
all electrons, bound and free, contribute to the total cross section. Daughter products from 
the absorption processes were not followed (e.g. positron annihilation photons from pair 
production and K-shell fluorescence photons from photoelectric absorption.) In order to 
test the validity of this last approximation, we compared the results of our 3-dimensional 
Monte Carlo transport code (Maverick) with the 1-dimensional code FASTGAM (Pinto 
& Woosley 1988a). Good agreement was found between the model spectra whether the 
daughter products of absorption were included in FASTGAM or not. This justified the 
decision to disregard the fluorescence and annihilation photons in Maverick. 

Figure 8, shows comparison spectra with the 1-dimensional FASTGAM code for 
model lOHMM (Pinto & Woosley 1988b). This is a logarithmic plot of photon flux 
(photons/second/MeV/cm^, assuming a distance of 60 kpc) across the energy range 0.3 keV 
- 4 MeV. The input model (lOHMM) is an artificially mixed version of the Woosley 1988 
lOH model progenitor for SN 1987A. The first comparisons between FASTGAM and 
Maverick (left panel of Figure 8) showed significant differences in the location of the hard 
X-ray fall off. This was due to an invalid assumption in FASTGAM of constant electron 
fraction for all species which contribute to photoelectric and pair production opacities. 
This assumption is not valid for hydrogen, and results in a lower calculated opacity for 
the lower energy X-rays. The right panel of Figure 8 shows the hard X- and 7-ray spectra 
calculated from the two different codes after the absorption opacity correction was made. 
There is good agreement across the spectrum to within the uncertainties of the Monte Carlo 
calculation. 



3.2. Hard X-ray and 7-ray Spectrum 

Figure 9 is a logarithmic plot of photon flux in units of photons/second/ Me V/cm^ 
across the energy range investigated with these simulations (0.3 keV - 4 MeV). We 
have placed this object at the distance of the Large Magellanic Cloud (60 kpc) for easy 
comparison with fiux data from SN 1987A observations. The 5 panels are spectra from the 
different time slices; in each panel, we plot the spectrum for the Symmetric model, along 
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with polar and equatorial views of the Jet2 model. It can be seen immediately that the 
hard X-rays emerge earlier from the ejecta with a global explosion asymmetry (Jet2 model). 
This holds regardless of viewing angle (pole versus equator) towards the explosion. 

The fact that the hard X-ray flux in the aspherical explosion model is larger than the 
symmetric explosion, regardless of line of sight, can be understood in principle from optical 
depth arguments. In Figure 10, we show a contour plot of density (outer contour) and 
^^Co number density (inner contour) for the Jet2 and Symmetric models at t = 150 days. 
Decay of ^^Co is the major source of 7-ray photons, so the inner contour essentially traces 
the surface of the emission region. The horizontal and vertical lines in Figure 10 represent 
lines of sight from the ejecta surface to the emission source and are labeled with the optical 
depth along that line of sight. The dominant opacity for the hard X- and 7-rays is Compton 
scattering off electrons and, since the density contours remain roughly spherical in both 
models, the optical depth from a given point to the ejecta surface is roughly constant. In 
the Jet2 model, the ^^Ni was mixed out to larger radii in the polar direction, so it makes 
perfect sense that we see enhanced emission over the Symmetric model spectrum for that 
viewing angle (the optical depth that the high energy photons must pass through is 10 in 
the symmetric model versus ~ 6 — 7 along the polar line-of-sight in the Jet2 model) . 

At a first glance, one might expect that, in the Jet2 explosion along the equatorial 
line-of-sight, the total flux should also be low (the optical depth from the nickel in the 
equator is also roughly 10). However, this material does not dominate the high energy 
emission seen along the equatorial line-of-sight. The material ejected along the poles has 
been mixed far enough out in the ejecta that the optical depth these high energy photons 
must travel through, even along the equatorial line-of-sight, is quite low (~ 7). It is this 
nickel which dominates the hard X-ray emission at all viewing angles. In fact, the optical 
depth from the "ends" of the ^^Co distribution does not differ very much between the polar 
view and the equator view (r = 6 & 7 respectively), which explains why the overall hard 
X-ray flux depends only mildly on viewing angle. 

For the later time slices, this mismatch in escaping emission from the "ends" versus the 
equatorial plane ejecta is less pronounced, and the equator view spectrum has comparable 
contributions from both regions. 



Table 2. High Energy Luminosities^ 



Model'^ Timr 3-31) ke\"' kr-V kc-V 500-11)1)1) ke\'<' 1000-381)1) kc-V 847 kcV Line'^ 1238 kr\' Linr'' 



Syn 



160 
200 
250 
300 
366 

150 
200 
260 
300 
366 

150 
200 
250 
300 
365 



3.716(0.138) X 10^ 
6.172(0.044) xlO^ 
1.844(0.006) X 10^ 
3.714(0.007) X 10^ 
4.629(0.010) X 10^ 

1.503(0.585) X 10^ 
8.262(1.177) XIO-'' 
2.783(0.164) X 10^ 
4.260(0.279) xlO^ 



2.008(0.162)xl0^'' 
3. 183(0.055) X 10^'' 
1.272(0.009) X 10^* 
2.996(0.011) X 10^* 
4.910(0.018) X 10^^ 

9.891(4.221) x 10^'' 
5.565(1.103) X 10^''' 
2.290(0.249) X 10^* 
3.996(0.460) X 10^* 



4.676(0.218) X 10-^ ' 5.667(0.435) X lO"" 



2.136(0.710) xlO^** 
1.160(0.139) xlO^'' 
3.461(0.181) xlO^'^ 
4.556(0.286) X 10'" 
4.151(0.208) X 10'" 



1.447(0.516) X 10^'^ 
8.387(1.527) xlO-''^ 
3.191(0.301) X lO-'"^ 
5.015(0.520) X 10-''** 
6.217(0.464) X 10-''** 



6.637(0.930) X 10^^ 
1.068(0.032) X 10^^ 
4.302(0.062) X 10^® 
1.043(0.007) X 10^^ 
1.866(0.011) X 10^^ 

2.812(1.573) XIO^'^ 
1.779(0.543) xlO^^ 
7.969(1.379) X 10^® 
1.476(0.262) X 10^^ 
2.286(0.268) X 10^^ 

4.529(2.289) XIO^'^ 
2.887(0.783) X 10^** 
1.186(0.181) xlO'''^ 
1.995(0.319) X 10-''° 
2.746(0.300) X 10^^ 



9.281(2.135)xl0^'' 
1.440(0.075) X 10^^ 
5.727(0.119) X 10^* 
1.367(0.015) X 10^^ 
2.446(0.025) X 10^^ 

3.877(2.897) x lO"''^ 
2.066(1.004) X 10^** 
8.943(2.822) X 10^* 
1.774(0.553) xlO^^ 



3.973(1.448) XlO'"^ 
6.686(0.498) X 10^® 
2.016(0.074) X 10^® 
4.267(0.086) X 10^^ 
6.660(0.131) X 10^^ 

1.200(1.089) X 10^** 
6. 171(4.649) xlO^^ 
2.700(1.366) X 10^^ 
5.049(2.500) X 10^" 



2.933(0.613) X lO''" 7.678(2.686) X 10' 
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7.679(4.313) xlO^''' 
4.747(2.002) xlO^* 
1.824(0.441) xlO^'' 
2.987(0.782) X 10^" 
4.048(0.722) X 10^® 



3.592(2.866) xlO^'* 
1.904(1.098) X 10^° 
5.588(2.098) X 10^^ 
8.274(3.336) X 10^^ 
1.010(0.312) X 10*° 



6.828(1.593) X 10*1 
1. 156(0.061) X 10*^ 
5.008(0. 103) X 10*^ 
1.362(0.014) XIO" 
3.016(0.027) XIO" 

2.748(2.547) X 10*2 
2.150(1.188) X 10*-'' 
8.641(2.862)xl0*^ 
1. 937(0. 612)xl0** 
3.864(0.674) X 10** 

5.486(3.883)xl0*2 
4.442(1.973) x 10*-' 
2.244(0.439) x 10** 
4.357(0.884) x 10** 
7.200(0.899) X 10** 



4.687(1.588) X 10*1 
9.622(0.632) X 10*2 
4.476(0.106) X 10*^ 
1.280(0.014) X 10** 
2.886(0.026) X 10** 

3.041(2.813) X 10*2 
1.244(0.942) X 10*^ 
7.455(2.908) X 10*^ 
1.712(0.572) X 10** 
3.624(0.634) X 10** 

3.041(2.820) X 10*2 
4.628(2.118) X 10*^ 
2.334(0.448) X 10** 
4.210(0.821) X 10** 
6.383(0.832) X 10** 



^Numbers in parentheses represent 1 — cr Monte Carlo uncertainties. 

^The input models include the Symmetric explosion model (Sym), the Jet2 explosion model for the polar (Jet2-P) and equator (Jet2-E) viewing angles 
'^Units are erg s~l 
^^Units are 7 s~l 
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3.3. ^-ray Line Profiles 

Although the overall hard X-ray emission shows little variation between pole and 
equator views, a detailed look at the 7-ray line profile shapes and strengths, for the 1.238 
and 0.847 MeV ^^Co lines, does reveal trends with viewing angle. Figures 11 and 12 show 
line profiles of these two ^^Co lines for both the Symmetric and Jet2 explosion models. 
The broadening of the line is caused by Doppler velocity shifts resulting from the spatial 
distribution of radioactive nickel in the homologously expanding ejecta. The 4 panels shown 
are for days 200, 250, 300, and 365 after explosion. The three lines in the Jet2 spectra 
represent different viewing angles through the ejecta (along the pole, the equator and an 
intermediate angle ~ 45°.) For the Symmetric spectra, we have plotted these same viewing 
angles. 

3.3.1. Global Asymmetry 

Both explosion scenarios (Symmetric and Jet2) show blue-shifted line profiles and this 
effect is most enhanced in the Jet2 model. These differences can be best understood by 
examining the physical effects which dictate the formation of the line profile edges. The 
blue edge to the lines is set by the maximum observed line-of-sight velocity of the ^^Co in 
the ejecta. Because the nickel was mixed further out (and achieves higher velocities) for the 
asymmetric explosions, the 7-ray line profiles viewed along the line-of-sight of the jet are 
much more blueshifted for the Jet2 model than the Symmetric model 

Note, however, that the line profiles in the asymmetric explosion model (Jet2) depend 
upon viewing angle. 90° off the jet axis, the nickel producing the line profiles arises from 
both fast-moving jet ejecta and the ejecta along the equator. Because the fast-moving 
ejecta is moving perpendicular to the line of sight, the 7-rays from this material are not 
blueshifted significantly, and the slow moving equatorial material also produces only a 
modest blueshift. Hence, the line profiles viewed along the equatorial line-of-sight are 
blueshifted much less than those viewed along the polar line-of-sight. 

The red edge of the fines is determined by the escaping emission from ^^Co with the 
smallest line-of-sight velocity in the ejecta. In a Symmetric model, this should be an 
indication of how deep into the ejecta we can see along a given viewing angle. However, 
there is a more pronounced effect at play in the asymmetric explosion models. As we 
mentioned in the previous section, much of the 7-ray emission for the equatorial view arises 
from the "ends" of the ^^Co distribution. This material has a very low line-of-sight velocity 
for an equatorial observer, since it is being ejected predominantly in the polar direction. 
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This allows for a significantly lower velocity red edge of the equator view lines, even though 
the optical depth profiles do not vary much between polar and equator viewing angles. 

Also note that the fine centroids become less blue-shifted with time. As the supernova 
expands, the opacity in the ejecta drops and emission from material located deeper into 
the ejecta (and thus at smaller radial velocities) becomes visible. This results in an overall 
broadening of the line, as well as a redward shift of the line centroid. 



This basic understanding of the fine profile edges addresses only global features of the 
emitted decay lines. To understand the range of possibilities that the line shape may take 
between those edges, we must concentrate on the spatial distribution of the radioactive 
elements in the supernova ejecta. 

Because the expansion is basically homologous after 100 days, the hne-of-sight velocity 
of a fiuid element in the ejecta is proportional to its distance above the mid-plane of 
the explosion. Each spectral energy bin in the line profile can be mapped to a unique 
line-of-sight velocity in the ejecta, which can in turn be mapped to a specific height above 
the mid-plane. So, given the distribution of cobalt as a function of x, y and z-coordinates 
( ■^Co{,x,y,z) , where z is along the line-of-sight to the observer), the line profile shape 
should be proportional to 



Optical depth effects, as discussed above, will alter the line shape from this basic 
picture, but it is important to realize that the clumping, which occurs as the inner material 
is mixed outward, can cause sizeable irregularities in the profiles. Figure 13 shows a plot of 
Nco versus line-of-sight velocity for the polar view of model Jct2 at 250 days with higher 
spatial resolution than the previous run. The emergent 7-ray line flux is plotted with a 
dotted fine and can be seen to peak at velocities of around 1800 km/s. This corresponds 
to the location in the ejecta at which the exponential attenuation from the optical depth 
surpasses the flux enhancement from increasing cobalt mass: tmcV ~ 4. As the ejecta 
expand and optical depth drops, the 7-ray line profiles will take on the irregular shape seen 
in the cobalt distribution. 



3.3.2. Clumping Asymmetry 




(3) 
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4. Conclusions 

In our simulations, we found that asymmetric explosions lead to extensive mixing of 
the supernova ejecta, placing the products of explosive burning well into the helium layer of 
the star. Even mild 2:1 asymmetries can mix 10% of the nickel out to the inner edge of the 
hydrogen layer. If such mixing occurred in weak explosions, as well as the strong explosions 
presented in this paper, these mildly asymmetric supernova explosions could explain the 
extensive mixing required in population III stars (Umeda & Nomoto 2002) and black hole 
binary systems such as Nova Scorpii (Podsiadlowski et al. 2002). 

Prom the Monte Carlo transport simulations, we find that the high energy fluxes at 
early times (~150 days) for the asymmetric explosion model is roughly 4-5 times larger 
than for the corresponding symmetric explosion model. This suggests that emission from an 
explosion with mixing and global asymmetry would be observed roughly 50-100 days earlier 
than a symmetric explosion model with mixing alone. In addition, the line profiles for the 
asymmetric explosion vary with viewing angle and are brightest for the polar view. Along 
this viewing angle, the integrated line flux is enhanced by a factor of 2-3 over the symmetric 
explosion, and the lines peak at roughly 3-5 times brighter than the symmetric explosion 
line profiles. The integrated flux enhancement for the equator view of the Jet2 model is 
of order 30 % over the symmetric model. The line centroids observed for the asymmetric 
explosion are shifted relative to those from the symmetric model (more to the blue for the 
polar view and slightly more to the red for the equator view.) At later epochs in the SN 
explosion, the line proflles should reflect the irregular shapes seen in the underlying cobalt 
distribution, which result from the clumpy nature of the outward mixing of nickel in the 
ejecta. Table 2. summarizes the calculated luminosities of the different models for various 
continuum energy bands and the 1.238 and 0.847 MeV ^^Co lines. 

With the recent launch of the International Gamma-Ray Astrophysics Laboratory 
(INTEGRAL), it is particularly interesting to analyze our results in terms of the 
instrumental resolution and sensitivities this satellite can provide across these energy bands. 
It is important to keep in mind that the nickel mass synthesized in this 15 Mq model is 
roughly 3 times larger than the typically observed value (~0.08Mq) for core-collapse SN 
explosions. Since the high energy emission should scale linearly with nickel mass, the fluxes 
observed for this model will be roughly 3 times brighter than we would expect for a more 
typical core-collapse event. In fact, the high energy continuum for our Symmetric model 
explosion is of order 3 times brighter than was observed towards SN 1987A. Since 1987A's 
nickel yield was fairly typical, this comparison seems to validate the linear scaling of flux 
with nickel mass. 

Our models suggest that 7-ray line observations will be most useful for diagnosing the 
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departure of a core-collapse explosion from spherical symmetry. At energies around 1 MeV, 
INTEGRAL will have a spectral resolution of 2 keV and a narrow line sensitivity (3 — cr in 
10^ seconds) of ~5xl0~^ phot s~^ cm~^ (Hermsen & Winkler 2002). Our model lines are 
about 5 times broader than this resolution element, so the sensitivity for detecting them is 
worse by roughly ^/5. Using these specifications, INTEGRAL would be able to detect the 
^^Co lines from the polar view of the Jet2 model at a distance of ~650 kpc. 

We would like to thank P. Pinto, T. Evans, T. Urbatsch, A. Hcger, K. Nomoto, S. 
Wooslcy, J. Cannizzo and N. Gchrcls for their advice and encouragement on this paper. The 
large hydrodynamical simulations were run on the Avalon cluster at LANL. This work was 
performed under the auspices of the U.S. Department of Energy by Los Alamos National 
Laboratory under contract W-7405-ENG-36. Funding was also provided by a Feynman 
Fellowship at LANL, DOE SciDAC grant number DE-FC02-01ER41176 and Beomax, Inc. 
These simulations were run on the Avalon and Space-Simulator clusters at LANL. 
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Fig. 1. — Velocity distribution of nickel (solid), Oxygen (dotted), Helium (dashed), and 
Hydrogen (long dashed) in our 3-dimensional simulations. Comparing these distributions 
to the 2-dimensional simulations in Fig. 1 of Herant & Woosley (1994), we note that the 
distribution of elements is similar in both the 2- and 3-dimensional simulations. However, 
the stronger 3-dimensional explosion causes all of the ejecta to be moving slightly faster than 
that of the 2-dimensional simulation and it is difficult to compare mixing instabilities. 



Fig. 2. — Density contour (7xl0~^ gcm"'^) plot of the early stages of the convection, 4.3 
hours after the launch of the explosion. Notice that tendrils mixing out the material have 
already developed. It is this mixing that places nickel far beyond its initial distribution. 




Fig. 3. — 3- dimensional simulation of the symmetric explosion 1 year after the shock launch. 
The contours represent the cobalt distribution with a number density of 10~^. The colors 
denote the density distribution. Note that although the explosion is symmetric, Rayleigh- 
Taylor instabilities mix out the nickel. 
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Fig. 4. — Distribution of the nickel ejecta in mass, comparing jet explosions with a symmetric 
explosion (top panel) and equatorial explosions with a symmetric explosion (bottom panel). 
Note that as we increase the degree of asymmetry (Jet4 versus Jet2, Equator4 versus 
Equator2), the mixing increases dramatically, placing nickel well into the hydrogen envelope 
of the star. The dashed line (top panel) shows the extent of mixing if all of the nickel/cobalt 
decay energy is deposited into the nickel ejecta and produces almost as much mixing in a 
symmetric explosion as the Jet2 model. 
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Fig. 5. — 3- dimensional simulation of the Jet2 asymmetric explosion 1 year after shock 
launch. As in Fig. 3, the contours represent the cobalt distribution with a number density 
of 10~^. The colors denote the density distribution. The nickel is mixed out extensively in 
the polar direction where the explosion was strongest. However, the density distribution did 
not gain large asymmetries and remained fairly symmetric. 
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Fig. 6. — Distribution of the nickel ejecta versus velocity comparing jet explosions with 
a symmetric explosion (top panel) and equatorial explosions with a symmetric explosion 
(bottom panel). Note that as we increase the degree of asymmetry (Jet4 versus Jet2, 
Equator4 versus Equator2), the mixing increases dramatically, producing nickel velocities 
in excess of 3000 km s~^. The dashed line (top panel) shows the velocities achieved if all of 
the nickel/cobalt decay energy is deposited into the nickel ejecta and produces almost as 
much mixing in a symmetric explosion as the Jet2 model. 
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Fig. 7. — Fraction of nickel ejecta mixed out into the star for our set of simulations. Note 
that for mild asymmetries (Jet2) or if decay energy is included in a symmetric explosion 
(Symmetric+Decay), nearly 10% of the nickel mass is injected into the hydrogen envelope. If 
this amount of mixing occurs in weak explosions, an explosion that resulted in a 4.5 Mgblack 
hole remnant would still eject a moderate amount of nickel. 
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Fig. 8. — Total hard X- and 7-ray spectrum comparison between the ID Monte Carlo 
transport code FASTGAM and the 3D code used in this work (Maverick). The left 
panel shows the two calculated spectra before a correction to the absorptive opacities in 
FASTGAM. The right panel shows the comparison once this correction had been made. The 
agreement between the two code results is quite good. 



- 28 - 



- ■ ■ Jets (pole) 




-2 -1 

Log Energy (MeV) 



Fig. 9. — Total hard X- and 7-ray spectrum at 5 different times during the explosion 
(150,200,250,300,365 days) for symmetric (sohd lines) and aspherical (Jet2) explosions 
(dotted and dashed lines). The flux is determined by assuming the object is 60kpc from the 
observer. The dotted lines refer to an aspherical explosion where the jet is directed along the 
line-of-sight of the observer. The dashed lines refer to an explosion where the observer line- 
of-sight is directed 90° off of the jet axis, in the equatorial direction. Regardless of observer 
viewing angle, the aspherical explosion is ~2 brighter than the symmetric explosion. 
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Fig. 10. — Contour plots in the xz-plane of the Symmetric and Jet2 explosion models at t 
= 150 days. Inner contour is for ^^Co number density which traces the surface of the 7-ray 
emitting region. Outer contour is for the mass density which follows electron density and 
thus traces the dominant opacity source (Compton scattering). The lines represent lines-of- 
sight through the ejecta for which the optical depth from emission region to ejecta surface 
has been calculated. Regardless of viewing angle, the optical depth of the ^^Co ejected along 
the poles in the Jet2 explosion remains quite low. Hence, it is this material that dominates 
the observed emission for all viewing angles in the aspherical explosion. 
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Fig. 11. — Line profiles of tlie ^^Co 1.238 and 0.847 MeV lines for the Symmetric model at 4 
different times during the explosion (200, 250, 300, 365 days). 3 different viewing angles are 
shown: polar view (dotted lines), equatorial view (dashed lines) and an intermediate view 
angle of ~ 45° (dash-dot lines). The line profiles do not show significant variation with 
viewing angle (as would be expected for a symmetric explosion.) 
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Fig. 12.— Line profiles of tlie ^'^Co 1.238 and 0.847 MeV lines for the Jet2 model at 4 
different times during the explosion (200, 250, 300, 365 days). 3 different viewing angles are 
shown: polar view (dotted lines), equatorial view (dashed lines) and an intermediate view 
angle of ~ 45° (dash-dot lines). The flux axis is scaled by a factor of 4 over the Symmetric 
model profiles shown in Figure 11. Significant variations in the line profiles with viewing 
angle are apparent, and can be explained by considering the velocity distribution (and thus 
radial distribution in a homologous expansion) of the ejecta responsible for the observed 
emission. 
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Fig. 13. — Number of cobalt nuclei (Nco('2^)) versus ejecta velocity along the polar axis for 
the Jet2 model at 250 days with roughly 3 times more spatial resolution than the previous 
run. The 7-ray line flux for this polar view is plotted as the dotted line for comparison. 
Positive velocities correspond to blueshifted energies. The line flux at this epoch departs 
from the cobalt distribution at velocities of order 1800 km/s due to optical depth effects. 
As the supernova expands the line profile will take on the irregular shape of the underlying 
distribution. 



